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Although effective for two dimensional (2D) systems, some approximations may fail in describing 
the properties of one-dimensional (ID) models, which belong to a different universality class. In 
this paper, we analyze the adequacy of the Composite Operator Method (COM), which provides a 
good description of many features of 2D strongly correlated systems, in grasping the physics of ID 
models. To this purpose, the ID Hubbard model is studied within the framework of the COM by 
considering a two-pole approximation and a paramagnetic ground state. The local, thermodynamic 
and single-particle properties, the correlation functions and susceptibilities are calculated in the case 
of half filling and arbitrary filling. The results are compared with those obtained by the Bethe Ansatz 
(BA) as well as by other numerical and analytical techniques. The advantages and limitations of 
the method are analyzed in detail. 



I. INTRODUCTION 

The physics of interacting electrons confined to ID sys- 
tems is one of the most interesting fields of research in 
Condensed Matter Physics. The reasons are various. On 
the one hand, the physics of ID systems challenges the 
standard picture of interacting electrons in metals, which 
has the Fermi liquid (FL) theory as its basic cornerstone. 
As Tomonaga and Luttinger 7 7 showed, a strictly ID 
interacting electron system cannot be described by FL 
theory. In such a system charge and spin degrees of 
freedom merge into collective low-energy excitations that 
propagate with different velocities and the quasi-particle 
picture, essential to FL theory, breaks down. This new 
electronic state is called Luttinger liquid (LL). Signals 
of LL behavior can be sought in any physical realization 
of ID electronic systems. Synthetic organic metals like 
the Bechgaard salts are probably the best candidates 7 . 
These metals have crystal structures consisting of alter- 
nating layers of organic donor molecules like TMTTF 
and TMTSF, and inorganic anions such as PF§ , SbFg or 
Br. Stacking planar molecules yield an overlap of the 
molecular orbitals that is greatest along the stacks and 
weaker between them, thus producing quasi- ID conduc- 
tors. Recent optical measurements in the metallic state 
of various Bechgaard salts have shown consistency with 
LL behavior 7 . On the other hand, the Bechgaard salts 
are subject of intensive studies 7 7 because they have a 
rich phase diagram, with antiferromagnetic, spin-Peierls, 
spin-density wave and superconducting ground states. In 
particular, the superconducting state shows some simi- 
larities to that of high-T c cuprates (high anisotropic con- 
ductivity, large and anisotropic critical field 7 and short 
coherence length 7 ). Also, the interplay between antifer- 
romagnetic and superconducting ground states and the 
strong sensitivity of T c to non-magnetic impurities indi- 
cate an unconventional superconducting mechanism that 



still remains to be determined 7 7 7 . All these properties 
are undoubtedly a strong motivation for better under- 
standing the physics of interacting electrons in ID sys- 
tems. One of the most suitable Hamiltonian to consider 
for this purpose is the ID Hubbard model 7 . This Hamil- 
tonian is exactly integrable by means of the BA 1 . In 
this way, many properties are known exactly within the 
numerics needed in the case of arbitrary particle den- 
sity and/or finite temperature. Namely, many ground- 
state properties 7 7 (total energy, local magnetization, 
magnetic susceptibility, etc.), charge and spin excitation 
spectra 7 , and some thermodynamic properties 7 7 7 can 
be exactly computed. However, the BA does not pro- 
vide a complete framework for describing the physics 
of the ID Hubbard model since many properties, like 
the correlation and spectral functions, cannot be evalu- 
ated from the BA wave function except for some limit- 
ing cases 7 7 (infinite interaction, static case, half filling) . 
Therefore, to compute these quantities, which are among 
the most relevant ones for describing real materials and 
getting a complete overview, we must consider other ap- 
proaches. Bosonization techniques, conformal field the- 
ory and quantum transfer matrix (qtm) and string the- 
ory investigations are analytic methods often used for 
ID models. They permit to compute key quantities like 
correlation functions, scaling relations between their ex- 
ponents and the velocities of spin and charge collective 
modes, but also thermodynamic quantities like the spe- 
cific heat and the charge and spin susceptibilities 7 7 7 . 
However, these methods need lengthy and complex cal- 
culations. The numerical techniques 7 7 7 7 7 7 7 7 , in- 
stead, are limited by the small size of clusters and the 
impossibility of reaching very low temperatures. We are 
thus interested in analyzing the adequacy of a simpler an- 
alytical calculation scheme for describing the physics of 
correlated electrons in ID models. This method, called 
the Composite Operator Method (COM), is based on 
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the choice of an appropriate combination of standard 
fermionic field operators as basis for describing the ex- 
citations of the system. The properties of the com- 
posite fields are self-consistently determined through the 
equations of motion, and the parameters that arise (in- 
ternal parameters) are used to fix the representation 
where the dynamics is realized 7 . This procedure recov- 
ers symmetries that are usually badly violated by other 
approaches 7 and provides a good description of many 
features of strongly correlated systems; it is in excel- 
lent agreement with numerical simulations on the local 
and integrated quantities 7 7 7 7 and explains success- 
fully some anomalous thermodynamic 7 7 and magnetic 
behaviors 7 7 observed in high-T c cuprate superconduc- 
tors. Nevertheless, approximations adequate in higher 
dimensions can fail when applied to ID systems. The 
BA provides a useful test for any approximate method. 
In this paper we study the adequacy of COM to describe 
the physics of the ID Hubbard model. We evaluate the 
local, thermodynamic and single-particle properties, the 
correlation functions and the susceptibilities. Our re- 
sults are compared to the BA ones, whenever available. 
We also discuss the agreement with other analytical and 
numerical techniques, in particular, the Renormalization 
Group (RG) and the quantum Monte Carlo (qMC). The 
advantages and limitations of the COM are discussed. 
Some preliminary results have already been published in 
Refs. ? and ? ; the present work provides an exhaus- 
tive overview of the application of the COM to the ID 
Hubbard model. 

The paper is structured as follows. In Sec. Q, the 
framework of the COM for the ID Hubbard model is 
extensively described. The model is presented, the basis 
chosen, the solution given and many physical quantities 
addressed. In Sec. HI, the results for half filling and arbi- 



trary filling are analyzed separately. Special attention is 
devoted to the case of quarter filling since, together with 
the half-filled case, it is believed to be the scenario for the 
Bechgaard salts 7 . Finally, in Sec. IV, some conclusions 



are given. In the Appendix, the two-pole approximation 
scheme is reported in some detail. 



II. THE METHOD 



interaction. The hopping matrix is given by 



where unitary lattice constant and only nearest neighbors 
are considered. We have fixed the energy scale in such a 
way that ta = 0. Hereafter, any energy will be presented 
in units of t and we will consider ft = &b = 1. 



B. The Basis 

In the case of the Hubbard model, a natural choice 
for the operatorial basis is the Hubbard doublet ^(i) = 
(£ (*)' ^ (*))> where 



ft (j) = c t (j) [i _ n (j)] rf(i) = c f (i) n(i) 



(3) 



These operators describe the atomic transitions at the 
site i (i.e., the transitions n — <-> 1 and n = 1 «-* 2, 
respectively). We have 



4io) = k) etjw 



|0) 



T)t \a) = (-r +1 |U) 4{ k> =0 (4) 

IU> 

where a, a =T (1) or j (2) and {|0), \a), |U)} is the 
vectorial basis on the single site. 



C. The Green's function and the COM solution 

Considering a two-pole approximation 7 (see App. ^) 
and a paramagnetic ground state, the Fourier transform 
of the single-particle retarded thermal Green's function 
G(i,j) = (TZ |*(i), *t(j)}) may be written as 



A. The Model 

The ID Hubbard model is described by the following 
Hamiltonian: 

H = Y,[ti3-»Sii\ J(i)c(j)+Uj2nT(i) n L (i) (1) 



where ct (i) = ^cj (i) , cj is the creation electron op- 
erator at the site i in spinor notation, n a (i) = c\ (i) c a (i) 
is the charge density operator for the spin a, \x is the 
chemical potential introduced to control the particle den- 
sity (i.e., the filling) n and U is the intrasite Coulomb 



o-W (k) 



■iu-Ei (k) + ie 
where the spectral functions a^(k) are given by 

a " (fc)=/l1 AQ{k) 



(5) 



Mm _ _ I m ( fc ) 

2Q(k) 



w 2Q(fc) + (-Yg(k) 

<>22 '( k ) = J 22 7^7T-, 



4Q(fc) 



(6) 
(7) 
(8) 
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and Ei{k) = R(k) — (— YQ(k) are the energy spectra, 
with 



i2(fc) = ic/- M -2*cosfc--^- 



Q(k) 



1 



lg 2 (k) 



^m?(k) 
I11I22 



m(k) 



U 



g(k) = (l-n) 

ill-' 22 

m(k) = 2t [A + [p - I22) cos k] 



(9) 
(10) 

(11) 
(12) 



In = 1 — n/2 and I22 = n/2 are the diagonal elements 
of the normalization matrix (see App. |X|). As we can 
see from Eqs. (|i^)-([r2"|), the Green's function depends on 
the model parameters t and U, the external parameters n 
and T (temperature), and three internal parameters: the 
chemical potential A and p. The latter two parameters 
have the following expressions 

a = <n*K f «) - <?f w^)) (13) 

P=\K®nM ((c T (i)c i (i)) 0, c{(i)4(i)) (14) 

and they are related 7 to the difference between the hop- 
ping amplitudes within the two Hubbard subbands (A) 
and the intersite charge, spin and pair correlation func- 
tions (p). The superscript a indicates the projection on 
the first neighbor sites 



(15) 



71^(1) = (i) a ^ c(i) are the charge (/1 = 0) and spin 
(/i = 1,2,3) density operators, where — (1,1?), 
ct^ = (— I, a) and a are the Pauli matrices. The main 
effect of the internal parameters A and p on the bands 
Ei{k) is a uniform shift and a bandwidth renormalization, 
respectively. Depending on how these internal parame- 
ters are fixed 7 , very different results are obtained. In 
the COM they are determined by solving the following 
system of coupled self-consistent equations, 



n = 2(l - C 1X - 
A = Cfi - C 2 « 2 
C12 = 



C22) 



(16) 



The 



withCV = (* 7 (i)*t(i)) andC« 5 = (*?(t)*J(i) 

first two equations come from the existing relations be- 
tween the parameters n and A and the Green's function 
matrix elements, whereas the third equation has been 
chosen in order to satisfy the Pauli principle at the level of 
matrix elements 7 7 . This request, which fixes the proper 
representation of the Hilbert space, naturally implies the 
fulfillment of the constrains coming from the particle- hole 
sysmmetry 7 ; i.e., 



ju(2 — n) = U — n(n) 
A(2-n) = -A(ra) 
p(2 — n) = 1 — n + p{n) 



(17) 
(18) 
(19) 



Let us note that, at half filling, the third self-consistent 
equation in ( |l6| ) is identically satisfied and the p param- 
eter must be calculated by analytic continuation. In this 
case we have /i = £7/2, A = and the energy spectra 
(i = 1,2) have the following expression 



Ei(k) = -4tp cos k—i-y^fu 2 + l&t 2 {2p- if cos 2 k 

(20) 

D. The physical properties 

Within this calculation scheme the evaluation of the 
physical properties is straightforward once the internal 
parameters are determined. In the following we will de- 
scribe how the relevant quantities can be computed. 

1. The local quantities 

The chemical potential /1 is one of the outputs of 
Eqs. ([l6j). In the non-interacting case, \i is detetmined 
as a function of the particle density n and of the temper- 
ature T by means of the equation 



n = 1 



1 



{l-2f F [E(k)]}dk 



where 

E(k) = —[i — 2t cos k 



(21) 



(22) 



is the non-interacting energy spectrum and /f(^) is the 
Fermi function. At zero temperature, the previous equa- 
tion can be solved analitically and gives 



-2t cos ( — n 



—1 

\2 



(23) 



The internal energy per site E can be calculated as the 
thermal average of the Hamiltonian and is given by 



7<5 



(24) 



where D — (n|(i) nj,(i)) is the double occupancy. In the 
insulating phase, at zero temperature and half filling, the 
previous equation assumes the following expression 



E= U VU 2 + at 2 ( 
4 7r(2p-l) 1\ 

_ U 2 (2p + 1) 

2tt (2p - 1) VU 2 + at 



at 2 



U 2 + at 2 



at 2 



U 2 + at 2 



(25) 



where a = 16 (2p — 1) , K(x) and £ {x) are the complete 
elliptic integrals of first and second kind, respectively. In 
the non-interacting case, we have 



1 r 

E = -4t~ / F [E(k)] cos k alk 
n Jo 



(26) 
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We recall that, at half filling, the Bethe Ansatz result of 
Lieb and Wu ? reads as 



E = -4t 



J (x)J 1 (x) 
o x(l + e xU / 2t ) 



dx 



(27) 



where J n (x) is the Bessel function of order n. The double 
occupancy D can be obtained as 



D = - - C 22 



(28) 



In the non- interacting case, we have D = The local 
magnetization Lq = j (113(1) n^(i)) can be computed by 

L ° = Uh D ) (2S » 

In the non-interacting case, we have Lq = |n (l — 

2. The thermodynamics 

The thermodynamic properties can be computed 
through the appropriate integrals of the chemical poten- 
tial and its derivatives. In particular, the specific heat C 
can be obtained as 



C(U,T,n) = -r£ ^n(JJ,T,ri)dn 



(30) 



where we have made use of the thermodynamic relation 
E = F+TS, with F and S the Helmholtz free energy and 
the entropy per site, respectively, and we have exploited 
the following expressions 



(31) 
(32) 



F(U,T,n) = / n(U,T,n')dn' 
Jo 

S(U,T,n) = -J —fi(U,T,n')dn' 
which give 



E(U,T,n) = 



^U,T,n')-T—^U,T,n') 



dn! 
(33) 



As it was shown for the 2D case 7 , the temperature 
derivatives of the chemical potential can be expressed in 
terms of the internal parameters. It remains clear that, 
once the self-consistent equations (|l6|) are solved, there 
exist two ways to calculate the physical quantities. On 
the one hand, we can exploit, whenever is possible, the 
relations with the Green's function matrix elements; on 
the other hand, we can use the relations between the con- 
jugate variables and the Helmholtz free energy computed 
through the chemical potential. The energy per site is a 
clear example of these two ways of calculation, since it 
can be computed by means of both Eqs. (p3) and ([33|). 



Another example is the double occupancy, which can be 
calculated by using Eq. ( p8| ) and as 



D(U,T,n) = -^F(U,T,n) = j 



d_ 

dU 



fi(U, T, n') dn' 



(34) 



Obviously, the exact solution of the model gives identical 
results whichever way we choose. On the contrary, any 
analytical approximation can receive different results be- 
cause the first way mainly exploits the computation of 
two-particle static correlation functions while the second 
way is based on the value of the chemical potential, which 
can be computed by using only one-particle static corre- 
lation functions; in this case, for particles, we intend the 
original c electrons. Another procedure to compute the 
Helmholtz free energy F and the entropy S exploits the 
relations E = F + TS and S = — Using the latter, 
we can rewrite the former as follows 



E _ d (F 
T^~dT\T 



(35) 



and obtain for <C 1 
F(n, T, U) = E(n,T*,U) 



T f T E(n,T*,U)-E(n,T',U) ^ 



T 



(36) 



where the value of the internal energy E is given by 
Eq. (^i|). Hereafter, we will put a subindex H near 
any quantity calculated from the matrix elements of the 
Green's function and a subindex T near any quantity 
computed through the value of the chemical potential. 



3. The single-particle properties 

The momentum distribution function n(k) is defined 
by means of the equation 



1 r 

= — / n(k) dk 



(37) 



Within the COM, it can be computed as follows 

n(k) = 2{f F [E 1 (k)]Z 1 (k) + f F [E 2 (k)]Z 2 (k)} (38) 
where 



Zi(k) = \ 



l-(-) 



(1 - n)g(k) 



2Q(k) 



(39) 



are the weights of the two subbands. It is easy to verify 
that Zi(k) < 1. In the non- interacting case, we have 



n(k) = 2/ F [E(k)\ 



(40) 



and, at zero temperature, the Fermi momentum fcp (de- 
fined by E(kp) — 0) assumes the value ^n. 
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The density of states for the c electrons is given by the 
following expression 



2 



In the non-interacting case, we have 



1 



1 



2vri 



1 - 



co + n 



2t 



(41) 



(42) 



which clearly exhibits the well-known ID van Hove sin- 
gularities at the edges of the band (cfr. Eq. ^2j). In the 
interacting case, each singularity splits in two as we have 
two distinct subbands. 



4- The correlation functions and the susceptibilities 

The distribution function B(r) = ^c T (r)c(0)) can be 
computed as follows 



B(r)= 28 rfi - 2^2 C ab {r) 



(43) 



where C a b{r) = I ^ a {r)W b {Qyj is the static correlation 
function given by 

C ab (r) = f E T e ikr {l-f F [E m (k)}}<r^(k)dk 

71 m =l "'-^ 

(44) 

At zero temperature and half filling, Eq. ([l3|) assumes 
the following simple expression 



s r . r~ - % 1 C cos (k r) cos k „ 

B(r = <*r,o - it (2 P - 1 — / dfc 

2vr J r(fc) 



(45) 



where T (k) = yU 2 /A + At 2 (2p - 1) cos 2 k. In the non- 
interacting case, we have 

i r 

B(r) = 5 r o - - / {1 - 2/ F [E(k)}} cos fcr dfc (46) 

T JO 



which, at zero temperature, becomes 

2 7r n r 

Sir) = — sin 

w Trr 2 



(47) 



showing damped oscillations (r x ) of wavelength A = 

4 _ 2tt 
n &f ' 

By considering the one- loop approximation' , the two- 
particle Green's functions can be calculated in terms of 
the single-particle ones. 



The density-density correlation function N(r) — 
(n(r) n(0)) will be computed as follows 

^)=« 2 + ^^r E ^Qa bac (r) (48) 



n-2D 



,.b,c=l 



where 

Qabac(r) = [Iab(r) - C ab (r)} C ac (r) 



(49) 



and I ab are the elements of the normalization matrix. 
It is worth noting that Eq. (|4|) satisfies the sum rule 
N(0) = n + 2D; this is a clear manifestation of the con- 
serving nature, at level of local sum rules, of the one-loop 
approximation. In the non-interacting case, the density- 
density correlation function reads as 



N(r) = n 2 + S r . n-^B 2 (r) 



(50) 



and shows damped oscillations (r~ 2 ) of wavelength A = 
— = at zero temperature. 

The spin-spin correlation function S(r) — (713(7") 713(0)) 
can be obtained from 



77 + 2D - 77 2 



(51) 



,b,c=l 



We can now establish the following relation between the 
spin and charge correlation functions 



S(r) 



n-2D 
n + 2D-n 2 



[N(r) 



(52) 



which implies they have the same spatial dependence. 
Consequently, within the one-loop approximation, gaps 
in charge and/or spin sectors will open simultaneously. 
This is also the situation in the non-interacting case, 
where the spin-spin correlation function reads as 



S(r) = 5 r , n-~B 2 (r) 



(53) 



It is worth pointing out that this limitation is connected 
to the use of the one-loop approximation and not to the 
Composite Operator Method. 

In the linear response theory, the spin susceptibility 
is defined as Xs{k,co)= — T (TZ [713(7-, t) 713(0)]) {T is the 
Fourier transform operator) and can be computed by 
means of the following expression within the one-loop 
approximation 



Xs(k,co) = 
where 



77(2-n) 

~n + 2D-n 2 



E laaQabacik^) (54) 



2 

i- £ { dpa { :i\k+p)cr£\p)x 

m.n—1 n 
fc[E m {p)\- fc[En{k+p)] 



lj + E n (k+p) - E m (p) + ie 



(55) 
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The spin susceptibility, in the non-interacting case, is 
obtained from the following integral 

ri \ If f F [E(p)}~f F [E(k+p)} 

Xs{k,uj) = / — — r— — dp (56) 

7T J_„ oj + E(p) - E(k +p) +ie 

At zero temperature, the uniform and static spin suscep- 
tibility reads as 



1 



Xs 



2N CC (0) 



(57) 



The charge susceptibility is defined as Xc(k,uj) 
T (7Z [n(r, t) n(Q)]) and it can be computed as follows 
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Once again, the spin and charge susceptibilities satisfy 
the following relation 



Xs(k,u) 



2D 



2D 



;Xc{k,uj) 



(59) 



In the non-interacting case, the charge susceptibility co- 
incides with the spin one (i.e., Xc(k,u>) = Xs{k,uj)). 

The uniform and static spin and charge susceptibilities 
can be also calculated directly from their thermodynam- 
ical definitions 



Xs 



dm 
~dh 
dn 

Xc TJ 



h=0 



(60) 
(61) 



where m and h are the magnetization per site and the 
external applied field, respectively. 

It is worth mentioning that all the expressions given in 
this Section, computed within the Composite Operator 
Method, exactly reproduces the non-interacting and the 
atomic limits. 



III. THE RESULTS 

The set of self-consistent equations in (|l6|), that deter- 
mines the internal parameters, admits two distinct solu- 
tions. Hereafter, we will call these solutions COM 1 and 
COM 2. The evolution of the internal parameters with 
the external ones (i.e., filling, intrasite Coulomb potential 
and temperature) reveals substantial differences between 
the two solutions. The A parameter is smaller and much 
more sensitive to the strength of Coulomb interaction in 
COM 1 (see Fig. |l|). For n < 1, the p parameter is neg- 
ative or very small and positive in COM 1, while it is 
always positive and of the order of the filling in COM 2 
see Fig. ||; actually, at half filling and on increasing U, p 
tends to in COM 1 and to 1 in COM 2. In COM 1, 
the chemical potential fi shows a discontinuity at half fill- 
ing for any finite value of the Coulomb interaction and 



FIG. 1: Internal parameter A as function of n for T = 0.01 
and U = 4, 8 and 12. 




1.75 2.00 

n 



FIG. 2: Internal parameter p as function of n for T = 0.01 
and U = 4, 8 and 12. 



a zone of instability (i.e., a negative compressibility) at 
small doping, temperatures and interaction strength (see 
Fig. ||). In COM 2, the discontinuity of the chemical po- 
tential appears only after a critical value of the Coulomb 
interaction is reached (see Fig. |]). As we will show in the 
next section, the absence of the Mott-Hubbard transi- 
tion, which is a consequence of the mainly negative value 
of the p parameter, plays a key role in the physics de- 
scribed by the COM 1 solution. The instability of COM 
1 can be studied by looking at the compressibility (i.e., 
k = ;j3-§^)- The instability is confined to a very small 
region in the plane (n,T). This region does not compre- 
hend half filling (see Fig. ||). For the 2D Hubbard model, 
which also has two solutions within the framework of the 
COM, the region of instability is much larger. 



A. Properties at half filling 

The ground-state properties of the ID Hubbard model 
at half filling were exactly derived by Lieb and Wu using 
the BA ? . This exact solution corresponds to an insulat- 
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FIG. 3: Chemical potential fj, as function of n for T = 0.01 FIG. 5: Instability region within the plane n-T for U = 4 and 
and U = 4, 8 and 12 (COM 1 solution). 8 (COM 1 solution). 
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FIG. 4: Chemical potential /i as function of n for T 
and U = 4, 8 and 12 (COM 2 solution). 
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ing state with short-range antiferromagnetic (AF ) corre- 
lations for any finite value of U. According to this, we 
will compare the evolution of the gap in the excitation 
spectrum and the band structure of both solutions with 
the BA exact results, in order to choose the solution that 
gives the most consistent physical picture at half filling. 

As already mentioned above, within the analysis of the 
chemical potential, the COM 1 solution presents a gap 
for any finite value of the Coulomb interaction, in agree- 
ment with the BA result, while the COM 2 solution is 
characterized by a critical value of the Coulomb inter- 
action (i.e., U c w 6.56) above which a gap opens (see 
Fig. H). The rate at which the gap opens in COM 1 
coincides with that of the exact solution for U > 4. 
The band structure (i.e., the excitation spectrum) of the 
two solutions is another interesting property to be stud- 
ied and compared. In fact, many features, more or less 
anomalous, of both solutions can be easily understood 
just looking at their spectra. As it can be seen in Fig. [t], 
the COM 1 solution has a typical AF band pattern 
(i.e., a quasi- halved Brillouin zone, the first excitation 
around k = and a very narrow bandwidth of the or- 



FIG. 6: Gap as function of U (COM 1 solution). 



der J = ^jf), in agreement with the BA result. On the 
contrary, the band structure of the COM 2 solution cor- 
responds to a typically paramagnetic state, with the first 
hole excitation at k — ±7r, the first electron excitation at 
k = and a bandwidth of the order 8t (see Fig. ||). In 
the figures, the energy is measured with respect to the 
chemical potential. In particular, while COM 2 has two 
subbands with both a minimum in k = and a max- 
imum in k = ±7r, the COM 1 solution has the upper 
subband with a maximum in k — ±7r and a minimum at 
k\ = ± arccos 9$ and the lower subband with a maximum 
at ki = 7r — k\ and a minimum in k = 0, where 



70 = 



pU 



2t(l - 2p)y/T=4p 



(62) 



For large values of the Coulomb interaction both k\ and 
ki tend to -| since p tends to zero (see Fig. ||). According 
to this analysis, the width of the subbands and the value 
of the gap in both solutions can be easily computed. In 
COM 2, the width of the subbands at half filling is W = 
8tp, which tends to 8t for large values of the Coulomb 
interaction, while the gap, above the critical value U c = 
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FIG. 7: Energy spectra E(k) for T = 0, n = 1 and U = 7 and 
10 (COM 1 solution). 
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FIG. 8: Energy spectra E(k) for T = 0, n = 1 and U = 7 and 
10 (COM 2 solution). 



4t^/4p — 1 w 6.56, has the expression 
A£ = -8tp+ v/f^ 2 + 16t 2 (2p- l) 2 



(63) 



On the contrary, in COM 1 the width of the subbands 
at half filling is 



1 



lr = 4/ 1, + ^VU 2 + 16t 2 (2p - If - U ^_ f (64) 



for large values of the Coulomb 



which tends to J 
interaction, while the gap, has the expression 



l-2p 



(65) 



Both expressions (i.e., Eqs. (|63j) and (pq)) for the gap 
tend to U for large values of the Coulomb interaction. It 
is worth pointing out that, within the two-pole approxi- 
mation, the p parameter rules the opening of the gap at 

half filling. In particular, if p < | ( 1 + we have a 

gapped solution, because the two subbands have oppo- 
site signs for any value of momentum and therefore they 



FIG. 9: ki and k%, see definition in the text, as function of U 
for T = and n = 1 (COM 1 solution). 




FIG. 10: Internal energy E as function of U for T 
n = 1 (COM 1 solution). 



and 



do not overlap. This is the case for the COM 1 solu- 
tion for any value of the Coulomb interaction. For the 
COM 2 solution we have a gap above the critical value 
U c . Otherwise, we have no gap and the two subbands 
overlap. In particular, both subbands have negative val- 
ues for \k\ < arccosxo, where xq = U/U c . 

We can conclude that only COM 1 gives a descrip- 
tion of the physics of the half-filled ID Hubbard model 
consistent with the exact results obtained by the BA. 
Therefore, in this section, we will mainly focus on this 
solution. 



1. Local properties 

The internal energy, at half filling and zero temper- 
ature, can be exactly calculated by means of the BA 
through Eq. (p7[). The two relevant limits (i.e., small 
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FIG. 11: Local magnetization Lq as function of U for T = 
and n — 1. 



and large Coulomb interaction) read as follows 



1*1 1 U 2 
]imE = -^ + ~U + 0(—) 
u^o it 4 t 



lim E = 



4-M2 



(66) 



where 4 In 2 = 2.77. The COM solution exactly agrees 
with the BA result in the weak-interacting limit, while 
in the strong-interacting limit, we have 



lim E 



COMl 
H 



t 2 r 



lim E§ OM2 



U 



U 3 ' 



(67) 



Again, while COM 1 gives a result very close to the BA 
one, the COM 2 solution is very far in this limit. The 
internal energy E at half filling and zero temperature, 
calculated by means of Eq. (|24|) in the COM 1 solution 
(i.e., Efl OM1 ), is shown, as a function of the Coulomb 
interaction strength, in Fig. [l(]. The results obtained by 
means of the BA 7 and other analytical approaches 7 7 7 
are also reported. As we can see, the agreement between 
COM 1 and BA is excellent. The self-consistent Lad- 
der approximation (SCLA) of Ref. ? shows also a very 
good agreement for all values of the coupling, but it does 
not have the correct behaviour for an infinite value of 
the Coulomb interaction. Moreover, both the Ladder 
{LA) 1 and the Gutzwiller (GWA and GWF) 1 approxi- 
mations go to zero at finite U , whereas the Renormaliza- 
tion Group (RG) ? has the right asymptotic behaviour 
for U i — *■ Co, but it does not reproduce the non- interacting 
limit. The local magnetization Lq at half filling and zero 
temperature as a function of the interaction strength is 
shown in Fig. nil We report both the COM results (i.e., 
Lqh M1 and L^? M1 ) and the results obtained by means 
of the BA 7 . L°° M1 is in excellent agreement with the 
exact solution. As one should expect, the electron local- 
ization increases with U and for infinite U reaches a sat- 
uration (i.e., zero double occupancy and zero kinetic en- 
ergy). Thus, the ID itinerant electron system described 
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FIG. 12: Specific heat Ct as function of T for n = 1 and 
U = 4, 8 and 12 {COM 1 solution). 
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FIG. 13: Local magnetization Lq as function of T for n = 1 
and U = 4, 8 and 12 {COM 1 solution). 



by the Hubbard chain is equivalent, at half filling and 
infinite U, to the system of localized spins described by 
the spin-i AF Heisenberg model. 



2. Thermodynamic properties 

The thermodynamic properties of the Hubbard chain 
can be evaluated using the BA by means of the finite- 
temperature formalism developed by Takahashi 7 . 

The specific heat Ct at half filling, calculated by means 
of Eq. Mia the COM 1 solution (i.e., Cf OM1 ), is shown 
in Fig. |12| as a function of the temperature for dif- 
ferent values of the Coulomb interaction. By increas- 



ing the interaction strength, the single peak present for 
U < W = 4t splits in two peaks moving to opposite di- 
rections in temperature (W is the non-interacting band- 
width). The low-T feature is due to spin excitations: it 
is located at T ~ J ; where J = is the magnitude 
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FIG. 14: Entropy S as function of T for n = 1 and U = 4, 8 

and 12 (COM 1 solution). FIG. 15: Static and uniform charge susceptibility Xc as func- 

tion of T (COM 1 solution). 



of the induced AF exchange parameter. Such a low-T 
peak is characteristic of the AF Heisenberg chain 7 . This 
feature is consistent with the AF-like band structure of 
the COM 1 solution (cfr. Fig. |). The high-T peak is 
obviously associated to charge excitations: it moves to- 
wards higher temperature as the gap increases with U. 
As T increases, the thermal energy allows the electrons 
to be excited across the gap. Such a structure of the 
specific heat, with low- and high- T regions, dominated 
by spin and charge excitations, respectively, is consis- 
tent with the physical picture described by the exact BA 
solution 7 7 7 7 7 and it also agrees with recent qMC cal- 
culations on finite chains 7 . 

In order to get information about the degree of local- 
ization of the electrons we have also computed the local 
magnetization as a function of the temperature. The re- 
sults are reported in Fig. [l^. Let us note that the max- 
imum localization occurs at finite temperature. This is 
due to the strong antiferromagnetic correlations present 
at zero temperature. These correlations require a virtual 
hopping and, therefore, diminish the degree of localiza- 
tion. Only by increasing the temperature we can sup- 
press the antiferromagnetic correlations and increase the 
localization. Note also that the flat regions at low tem- 
peratures end at the same temperature where the local 
minima develop in Ct- This further confirms the spin- 
nature of the low-T peak in the specific heat. 

In the temperature evolution of the entropy at half 
filling, spin and charge degrees of freedom manifest sep- 
arately in the strong-interacting regime U > W (see 
Fig. |l| where S% OM1 is calculated by means of Eq. B2| in 
the COM 1 solution). The low-T region is dominated by 
spin excitations and the high-T region by particle excita- 
tions. The internal energy shows an analogous behavior. 
The borderline between these two regions can be set at 
T ~ t in very good agreement with both BA' 1 and 
numerical 7 results. It is worth noting that, in our solu- 
tion, the entropy shows the correct limiting behavior for 
high temperatures (i.e., limT^oo St — 2 In 2 = 1.39). 





1 1 1 


1 1 1 1 






COM 1 


5 




n = 0.8 






U = . 


4 




U = 4 . 






U-8 






U = 12 


3 




2 
1 


-,\ 







i.i.i. 



0.0 0.5 1.0 1.5 2.0 

T 



FIG. 16: Static and uniform spin susceptibility \s as function 
of T ( COM 1 solution) . 



3. Correlation functions and susceptibilities 

A quantity which gives information about the evolu- 
tion of the charge excitations in the system is the charge 
susceptibility (xc), which we calculate by means of Eq. 
(|6l|). This property is shown in Fig. [l^ for n — 0.9 and 
ri=l; string theory results are taken from Rcf. 7 . The 
first thing to note is the strong dependence of Xc on the 
particle concentration at low temperatures and for any 
coupling regime (not shown): the charge susceptibility 
is strongly enhanced as n approaches half filling, while 
it goes to zero at n — 1 as a consequence of the open- 
ing of the gap in the charge excitation spectrum. It is 
worth noticing our good agreement with the practically 
exact results of qtm' . For increasing Coulomb interac- 
tion (not shown), the low-temperature enhancement of 
Xc in the low doping region is more evident, indicating 
that the charge excitations are strongly renormalized by 
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have obtained for the spin-spin correlation function (see 
Fig. |l7| ). It presents a typical paramagnetic behavior 
(i.e., it is always negative at any finite distance) in con- 
trast with the exact diagonalization results 7 7 which re- 
port short-range antiferromagnetic correlations (i.e., a 
spin-spin correlation function alternating in sign between 
sites and rapidly decreasing). Nevertheless, the present 
approach preserves the 2kp oscillations for any value of 
the Coulomb interaction. 



B. Properties at arbitrary filling 



FIG. 17: Spin correlation function S(r) for T = 0, n = 1 and 
U = 4, 10 (COM 1 solution). 

the Coulomb interaction near the metal-insulator transi- 
tion. At higher temperatures, however, \c decreases with 
increasing U regardless of the electron concentration 7 . 

These are already well-known results 7 7 that the ap- 
proximation considered here is able to reproduce, thus 
capturing the physics of the charge excitations near the 
metal-insulator transition. The temperature at which the 
gap closes due to the thermal excitations is somewhat 
larger in COM 1 than in the Bethe Ansatz results 7 7 . 
This is coherent with the larger value that we obtain for 
the gap. Also, \ c is more enhanced near half filling (for 
instance, at U = 8t and n = 0.8, Xc,mux ~ °- 9 while 
Xc,max ~ 0.25) because of the faster rate at which the 
gap opens in our approximation with respect to BA. 

Information on the physics of the spin excitations can 
be extracted from the evolution of the magnetic suscep- 
tibility with temperature and Coulomb interaction. We 
discuss the situation close to half filling. The magnetic 
susceptibility Xs is calculated by means of Eq. (|54j) and 
its behavior, as a function of the temperature and U, is 
shown in Fig. \s presents a peak at low tempera- 
tures which moves to lower ones as U increases. It is also 
worth noticing that \s is generally enhanced by U. At 
half filling [not shown], our results indicate a renormal- 
ization of the spin excitations, as it is obtained in the 
case of the charge susceptibility. This behaviour of Xs 
for the interacting half-filled chain does not reproduce 
the BA result 7 7 . According to the exact solution, Xs 
should show the same qualitative behaviour regardless of 
the particle density, because the metal-insulator transi- 
tion at n = 1 does not renormalize the spin excitations. 
Such disagreement between COM 1 and BA is due to the 
computation of the two-particle Green's functions within 
the one-loop approximation. In this approximation, the 
charge and spin correlation functions have the same spa- 
tial dependence (cfr. Eq. fl59[)), and consequently, gaps in 
charge or spin sectors will open simultaneously. There- 
fore, we obtain the same qualitative behaviour for both 
the charge and spin susceptibility. 

This disagreement is also present in the results we 



BA predicts a non-magnetic metallic ground state for 
the ID Hubbard model at arbitrary filling, with gapless 
charge and spin excitations 7 . Ground-state properties, 
like the internal energy, chemical potential and local mo- 
ment, were studied following Lieb and Wu as a function of 
the electron density and Coulomb interaction 7 7 , mainly 
in the large- U limit, where analytic expressions could 
be obtained. By means of the formulation developed 
by Takahashi for finite temperatures 7 , some thermal 
properties at arbitrary filling were also calculated" " . 
More recently, important single-particle properties, as 
the spectral function and momentum distribution func- 
tion, have been evaluated within the BA' 7 and by 
means of numerical techniques, like qMC 1 7 7 . Infor- 
mation on the charge and spin dynamics is the object of 
the most recent studies on one-dimensional models, with 
the aim of applying them to real systems. Thus, some 
authors have analyzed the corresponding two-body cor- 
relation functions of the ID Hubbard model by means 
of numerical approaches 7 7 , by using the exact BA 
solution 7 7 or through other analytic approaches, like 
g-ology 7 . 

In the following sections we present the results ob- 
tained for such quantities within the COM for the Hub- 
bard chain away from half filling. The COM 2 solution 
is not considered, since, as we have commented above, 
it does not provide a good description of the system in 
the case of half filling. Thus, the COM 1 results are 
compared to the ones available by other systematics. We 
devote special attention to the quarter-filled case, since 
it seems to be a relevant particle density for many ID 
organic metals 7 . 



1. Local properties 

The doping dependence of the internal energy, calcu- 
lated by means of Eq. (B4T) in the COM 1 solution (i.e., 
E c° M1 ) , is shown in Figs. |l8| and |l^ for two values of the 
Coulomb interaction. For comparison, we also report the 
BA results 7 and the ladder (LA) and self-consistent lad- 
der (SCLA) approaches 7 7 . COM 1 agrees reasonably 
with BA, reaching the best agreement at half filling. The 
ladder approximation 7 deviates more and more from the 
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FIG. 18: Internal energy E as a function of n for T = and 
U = 2. 



FIG. 20: Chemical potential (i as a function of n for T = 
and [7 = 2 and 4. 
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FIG. 19: Internal energy Basa function of n for T = and 
U = 4. 



FIG. 21: Double occupancy D as a function of U for T — 
and n = 0.5 and 0.75. 



BA as approaching half filling; the self-consistent ladder 
approximation 7 probes excellently at any doping. 

The evolution of the chemical potential with the parti- 
cle density for COM 1 is compared in Fig. ^o] with some 
numerical data 7 . There is a good agreement over all 
the range of doping for the higher temperature. The dis- 
agreement at the lower temperature can be understood 
by looking at the size of the gap in COM 1: COM 1 
gap is larger than the BA one and forces the chemical 
potential to assume lower (higher) values than in the BA 
solution for n < 1 (n > 1). 

The good agreement with BA, at quarter-filling and for 
any value of U, occurs for the double occupancy too. For 
other values of filling the agreement is less satisfactory 
(see Fig. H). 



2. Thermodynamic properties 

The evolution of the specific heat Cr with the particle 
density is studied for various temperatures in the strong 
coupling regime. As shown in Fig. E3, 7 = Cr/T has a 



peak at low densities that simply reflects the shape of the 
density of states (see discussion on single-particle prop- 
erties in the next section). As n approaches half filling, 
another peak develops at low temperatures. This indi- 
cates an increasing number of excited states due to the 
renormalization of the charge fluctuations at the opening 
of the Mott gap. Let us note that such a peak does not 
appear in strongly correlated systems which do not have 
a metal-insulator transition like the ID electron gas with 
delta-function interactions 7 . As one should expect, with 
increasing temperature, the two peaks merge. COM 1 
results for Ch /T recover qualitatively those obtained by 
BA 7 . There are some quantitative differences, namely: 
the first peak appears at higher fillings and the double- 
peak feature survives up to higher temperatures in BA; 
the peak near half filling is higher in COM 1. This latter 
difference is expected as the COM 1 gap is larger than 
the BA one. 

To complete the above discussion we show in Fig. ^ 
the particle density n versus the chemical potential /1 for 
various temperatures. COM 1 results are compared with 
the BA ones of Ref . ? . The agreement is very good at low 
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FIG. 24: Energy spectra _E(fc) for n = 0.5, T = and (7 = 2, 
2.5, 5 and 10 (COM 1 solution). 
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FIG. 23: Filling n as a function of the chemical potential fi 
for U = 8 and T = 0.05, 1, 2 and 4 (COM 1 solution). 



FIG. 25: Density of states for n = 0.5, T = and U = 2, 5 
and 10 (COM 1 solution). 



temperatures for densities smaller than 0.55. In the half- 
filled chain, T ~ t is a relevant temperature as it signs 
the border between T-regions dominated by either spin or 
charge correlations 7 7 . The agreement between COM 1 
and BA at T > t is very good for the whole range of 
filling. Of course, at higher temperatures COM 1 result 
reaches an excellent agreement with BA since the effect 
of correlations is completely suppressed. 

3. Single-particle properties 

The evolution of the band structure with the inter- 
action strength U is remarkable. Away from half fill- 
ing the AF-band pattern of the upper Hubbard subband 
in COM 1 disappears as U increases; the first electron 
excitations appear at k = ±7r and the first hole exci- 
tations move slightly away from the half-filling position 
k = ±tt/2 (see Fig. |J). 

The corresponding density of states is shown in Fig. [2^. 
It is well known that the density of states in the non- 



interacting case exhibits two van Hove singularities at 
the edges of the non-interacting band. In the interact- 
ing case each of these singularities splits in two since we 
have two distinct Hubbard subbands. Due to the ^4F-like 
band shape, three van Hove singularities appear in each 
subband leading to the six peaks of the figure. The third 
structure in the upper subband of Fig. ^5] is smoothed be- 
cause, away from half filling, as U increases the ^4-F-band 
shape of the excitation spectrum disappears. 

As we can deduce from the above results the system, 
away from half filling, is a conductor for any value of the 
interaction strength U in agreement with the BA results. 
Then, a natural question arises: what universality class 
does this conductor belong to? This is a central issue in 
the physics of ID '-models. 

The momentum distribution function n(k) is a rel- 
evant property because from its behavior the Fermi- 
liquid or non-Fermi liquid nature of the excitations can 
be inferred. As it is well known, a finite jump in the 
zero-temperature momentum distribution function at the 
Fermi momentum fcp would indicate that the quasiparti- 
cle excitations can be described by a Fermi liquid. Some 




FIG. 26: Momentum distribution function in real space B(r) 
for n = 0.5, T = and (7 = 0, 2.5, 5, 10 and 30 (COM 1 
solution) . 



variational 7 7 , analytic 7 and numerical 7 approaches 
failed in understanding the nature of the discontinuity. 
More recently, approaches valid in the weak coupling 
regime (g-ology and qMC) have found a power-law sin- 
gularity at fcp predicting a marginal Fermi-liquid nature 
away from half filling 7 . The power-law coefficient seems 
to be an increasing function of U and a decreasing func- 
tion of the particle density n. 

The approximation that we use limits the information 
that can be obtained from the momentum distribution 
function. Namely, any two-pole approximation has a mo- 
mentum distribution function with sharp discontinuities 
at the values of momentum where the Hubbard subbands 
cross the Fermi level. Despite this strong limitation, the 
momentum distribution function obtained at quarter fill- 
ing in COM 1 presents, besides a sharp jump at fcp, an- 
other discontinuity near 3/cf ■ This latter feature is also 
obtained, as a weak singularity, in BA calculations in the 
large- U limit 7 . 

The momentum distribution function in real space ob- 
tained in COM 1 is shown in Fig. ^6] for various in- 
teraction strengths and quarter filling. In the weak 
coupling regime, the distribution function has nodes at 
r = 4a, 8a, 12a, . . . which correspond to an oscillation of 
wavelength 8a. Since for these values of interaction and 
particle density the Fermi momentum is this is just 
a fcp oscillation. Such an oscillation is also observed by 
exact diagonalization results 7 for small U. For stronger 
interactions, a weak incommensurate modulation, in ad- 
dition to the kp oscillation, appears and corresponds to 
a 3&F oscillation again in agreement with the exact di- 
agonalization calculations 7 . The origin of this 3/cp fea- 
ture is however not clear. According to other numerical 
calculations 7 the 3fcp singularity seems to be a finite- 
size effect and vanishes in the thermodynamic limit. We 
actually confirm its presence also in the bulk system. 



FIG. 27: Spin correlation function S(r) for n = 0.5, T = 
and U = 2.5, 4, 5, 10 and 30 {COM 1 solution). 
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FIG. 28: Spin correlation function S(r) for U = 4, T = and 
n = 0.5 and 1 (COM 1 solution). 



4- Correlation functions and susceptibilities 

The spin correlation function gives information about 
the spin dynamics of the system. Combining this in- 
formation with the one obtained from the momentum 
distribution function, we can analyze how the electron 
dynamics is affected by the surrounding spins, since the 
spin configuration is modified when an electron moves to 
another site. Therefore, the electron motion will depend 
on how the spin configuration constrains the electron dy- 
namics. 

The spin correlation function in real space is evaluated 
in COM 1 for various interacting regimes and relevant 
particle densities. The results are shown in Figs. ^ and 
B8L As we can see, the amplitude of spin correlation in- 
creases with U and is smeared out away from half filling. 
The larger the Coulomb interaction is, the faster the cor- 
relations decay. The physical picture that emerges for 
the half-filled and quarter-filled Hubbard chain is differ- 
ent. At n = 1 there are strong antiferromagnetic spin 
correlations that decay very fast (at second neighbor sites 
they are almost zero). The existence of short-range AF 
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FIG. 30: Spin susceptibility Xs(k) for U = 4, T = 0.001 and 
n = 0.99 on the extended zone (COM 1 solution). 

order is also visible in the band spectrum as we previ- 
ously discussed in detail. For n = 0.5 these correlations 
are much weaker and decay more slowly than in the half- 
filled case. These results agree with the ones described in 
Ref. ? using exact diagonalization techniques. In these 
calculations 7 the spin correlation function has a 2kp os- 
cillation that is not smeared out away from half filling. 
Such oscillation is not observed in Figs. [27] and ^8] due 
to the very fast decay of the correlation amplitude, but 
it appears as a 2kp feature in the momentum-dependent 
spin correlation function S(k) (see Fig. [2^ and explana- 
tion below). 

The spin correlation function in k space has been re- 
cently studied by means of various approaches. For in- 
stance, by qMC in the weak interacting regime 7 7 , by 
means of BA in the strong interacting regime 7 and for 
any value of U by exact diagonalization technique 7 . All 
of them find a very narrow 2/cp peak which is incommen- 
surate away from half filling. These calculations were 
performed on finite-size systems. A detailed study of the 
size dependence of S(2kp) is given in Ref. ? , where it 
is shown how the peak narrows and increases with the 
system size. 



FIG. 32: Interband spin susceptibility xi n ' er (ft) for [7 = 4, 
T = 0.001 and n = 1, 0.99, 0.98, 0.96, 0.94, 0.92 and 0.9 
(COM 1 solution). 



The spin-correlation function in k space is shown in 
Fig. ^ for the COM 1 solution at quarter filling and var- 
ious interaction strengths. As U increases, an evolution 
towards a peaked curve at 2fcp is observed, in qualitative 
agreement with previous calculations 7 7 7 7 . We must, 
however, remark that within the present approach the 
peak at 2kp is much broader. At half filling, the height 
of the 2fcp peak is very much enhanced, in agreement 
with BA calculations in the large-?/ limit. The reduction 
of the S'(2/cf) peak as n decreases from half filling is due 
to the presence of holes moving in the system. 

The magnetic or spin susceptibility \ s gives informa- 
tion about the physics of the spin excitations. This prop- 
erty is calculated by means of Eq. (|4|) and analyzed as 
a function of temperature, momentum and interaction. 

Close to half filling we have a paramagnetic solution 
with period 2tt, but a strong AF order, with a quasi- 
halved Brillouin zone, is present (see Fig. ^0|). When we 
move away from half filling the central peak opens in two 
separate peaks (see Fig. [3l])). The incommesurability 
amplitude increases linearly with doping with a coeffi- 
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FIG. 33: Static and uniform spin susceptibility \s as function 
of T for n = 0.5 and U = 2, 3 and 4. 



FIG. 35: Spin susceptibility x B (k) for n — 0.5, T = and 
J7 = 0, 2, 4 and 10 (COM 1 solution). 





FIG. 34: Static and uniform spin susceptibility \a as function 
of T for U = 8 and n = 0.5, 0.7 and 0.9 (COM 1 solution). 



FIG. 36: Static and uniform spin susceptibility Xa as a func- 
tion of U for n = 0.5 and T = 1/14.5. 



cient of tt/2 at zero temperature. If we distinguish the 
intraband and interband contributions we can show that 
the interband contribution is very little (see Fig. |32|). At 
half filling the susceptibility is strongly reduced and goes 
to zero at zero temperature: the interband contribution 
goes much slower to zero with temperature and for the 
analyzed temperature (T = 0.001) is the only one present 
(see Figs. |3l] and p2| ). 

Fig. [33| shows the temperature evolution of \s for weak 
and intermediate coupling. Monte Carlo data 7 and qtm 
results 7 are included for comparison. The agreement is 
good, with a maximum deviation of 15%. We shall note 
that the position of the maximum does hardly depend on 
the strength of interaction and is located around 0.2 (in 
units of t). This indicates that the energy necessary to 
excite the spin modes does not depend on U. Since U is 
of Coulomb origin and affects the charge degrees of free- 
dom, the U -independence of the peak in the spin suscep- 
tibility shows that the spin excitations are independent 
of the charge. This is coherent with a Luttinger liquid 
description of the quarter-filled Hubbard chain. On the 



contrary, at half filling the position of the T-dependent 
spin susceptibility does strongly vary with U (see Fig. |l^ 
and Ref. ? ), indicating the breaking of the Luttinger 
liquid picture. 

Fig. 34 shows the strong coupling spin susceptibility 
for several particle densities. Once again, COM 1 results 
at quarter filling are in good agreement with BA 1 7 (see 
value and position of x" lax ); whereas it severely disagrees 
in the case of higher particle-densities. 

The low temperature, momentum dependence of the 
magnetic susceptibility Xs{k) at quarter filling is shown 
in Fig. |35| for all coupling regimes. As U increases, the 
peak moves towards lower k, indicating that spin excita- 
tions of larger wavelength mostly contribute. Also, two 
satellite structures appear and increase with U. Such 
structures have their origin in the van Hove singularities 
of the density of states. They are not observed in qMC 
studies probably because they exist only at very small 
temperatures (see Fig. |38| ). 

The static uniform susceptibility is plotted in Figs. 3q 
and E37| as function of U and T, respectively. In Fig. 36 
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FIG. 37: Static and uniform spin susceptibility Xs as a func- 
tion of T for n = 0.5 and (7 = 4. 
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FIG. 38: Normalized spin susceptibility Xs(k)-^ for n = 0.5, 
U = 2 and T = 0.005, 1/14.5, 1/7.25, 1/4.5, 1/3.75 and 1/3. 

we also report the quantum Monte Carlo data ? and the 
BA results. The agreement, in particular in the weak 
coupling regime, is quite good. In Fig. ^37] the compar- 
ison with the qMC 1 and the qtm 1 results shows: 1) 
the qMC data do not agree at all with the qtm results, 
which are practically exact; 2) a good qualitative agree- 
ment between COM and qtm at very low temperatures 
for all three values of Coulomb repulsion U( the qtm re- 
sults present a more pronounced peak) and a very good 
quantitative agreement at intermediate and high tem- 
peratures. Anyway, the position of the peak is very well 
reproduced confirming once more the capability of the 
present method to catch the spin energy scale although 
the retained spin correlations are weaker than what is 
expected according to the exact and numerical results. 
The weak coupling spin susceptibility (normalized to 

its non-interacting value at k = 0: Xs(0)[U = 0] = 
see Eq. ^) is shown in Fig. |38|. Despite the simplicity 
of the one-loop approximation we get a good qualitative 
agreement with the quantum Monte Carlo data 7 , as op- 



posed to RPA fits, which need to take different values for 
the renormalized interaction, as it is remarked in Ref. ? . 
We exactly reproduce the peak at 7r/2, which vanishes on 
increasing temperature, and the asymmetry they found 
in the intensity between k = and k = tt. 



IV. CONCLUSIONS 

In this paper we have analyzed the adequacy of the 
COM to describe the physics of correlated electrons in 
ID systems; in particular, we studied the ID Hubbard 
model. Various physical properties, like ground-state and 
single-particle properties, the thermodynamics, suscepti- 
bilities and some correlation functions, have been calcu- 
lated and compared to results obtained by means of the 
Bethe Ansatz, when available, and other analytic and 
numerical techniques (qMC, GA, SCLA, RG, g-oloqy). 

By considering a two-pole approximation and a param- 
agnetic state, the model is solved within the COM. Two 
mathematical solutions (COM 1 and COM 2) are ob- 
tained. In the case of half filling an analysis of the pa- 
rameters and the energy spectra that characterize the so- 
lution shows that only COM 1 is consistent with the BA 
results. Hence, the subsequent properties are discussed 
only for this solution. The half-filled and arbitrary-filled 
cases are addressed separately since they lead to different 
characteristics. 

The essential physics of the half-filled Hubbard chain 
is reproduced satisfactorily by the COM 1 solution. A 
gap opens for any finite value of the U interaction and 
the bands show the characteristic AF -like features, with 
bandwidth of the order of the antiferromagnetic exchange 
interaction J = This indicates, therefore, an insula- 
tor with short-range AF correlations, as it is known from 
the exact BA result. Also, the evolution of the total en- 
ergy, the double occupancy and the local magnetization 
with U are in excellent agreement with BA and improves 
significantly the results of other analytic approaches. 

The thermodynamic properties give the following pic- 
ture. As the interaction strength increases relative to the 
non-interacting bandwidth W and reaches U > W = At, 
two energy scales manifest in the system in the form of 
two peaks in the specific heat. These peaks are located 
at low (T ~ J) and high T and are associated to spin 
and charge excitations, respectively. The spin nature of 
the low-T peak in C't is confirmed by the evolution of the 
local magnetization with temperature. In the strong in- 
teracting regime, the spin and charge degrees of freedom 
also manifest separately in the entropy, where T ~ t is 
the border between the two regions. This picture agrees 
qualitatively and quantitatively (position and height of 
peaks in Ct and border between the spin and charge- 
dominated regions in St) with BA. 

The behaviour of the charge excitations near the metal- 
insulator transition is also very well captured. This can 
be seen from the analysis of the charge susceptibility Xc 
versus T for particle densities approaching n = 1. The 
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agreement with BA is only qualitative because the faster 
opening of the gap in COM 1 leads to larger values of 

Xc- 

The physics of the spin excitations, extracted from the 
evolution of the spin susceptibility Xsi is not properly 
described. Our results for Xs indicate, in disagreement 
with BA, a renormalization of the spin excitations as the 
system approaches half filling, analogous to what is ob- 
served for the charge excitations. This failure is inherent 
to the one-loop approximation; within the latter gaps in 
charge and spin sectors open simultaneously since the 
charge and spin correlation functions have the same spa- 
tial dependence. 

The basic physics of the arbitrary-filled Hubbard chain 
is reproduced satisfactorily by the COM 1 solution. In 
agreement with BA results, we obtain a non-magnetic 
metal for any coupling regime. Namely, the band spec- 
trum is gapless for any value of U, it loses gradually its 
AF-like characteristics, and the amplitude of the spin 
correlation function is much reduced in comparison with 
that obtained at half filling. 

Our analysis for arbitrary filling is centered mostly on 
n = 0.5 which is relevant to real quasi-i-D systems. At 
this particular filling, low-temperature local quantities 
like the double occupancy and the chemical potential are 
in very good agreement with BA results for both weak 
and strong coupling regimes. The internal energy versus 
particle density in COM 1 has a reasonable agreement 
with BA although it is not so good as that of other ap- 
proaches. 

The evolution of the specific heat with particle den- 
sity reproduces qualitatively the BA results. The In- 
dependent spin susceptibility at quarter filling has a 
reasonable agreement in the weak and strong coupling 
regime with quantum Monte Carlo and BA data, respec- 
tively. 

The momentum distribution function n(k), the corre- 
lation functions and the susceptibilities provide informa- 
tion on the universality class this system belongs to and 
on its charge and spin dynamics. 

We can grasp some characteristic features. Namely, 
the weak singularity at 3/cf obtained in the large-?/ limit 
BA calculations at quarter filling would correspond to the 
discontinuity near 3fcp that is observed in the COM 1 mo- 
mentum distribution function. This feature would also 
agree with the 3fcp oscillation of the distribution func- 
tion in real space C(r) obtained by numerical techniques. 
The kp oscillations of C(r) observed by exact diagonal- 
ization calculations are also reproduced in the weak cou- 
pling COM 1 results. 

The approach manages to grasp the different physics 
for the half-filled and arbitrary-filled case, in particular 
quarter filling. By comparing the amplitude of the spin 
correlation function in real space S(r), we conclude that 
at quarter filling these correlations are much weaker and 
decay more slowly than at half filling, in agreement with 
exact diagonalization results. 

The 2kp singularity of the fc-dependent spin correla- 



tion function S(k) obtained in recent BA and numerical 
approaches is qualitatively described in COM 1 at quar- 
ter filling, where, by increasing U, an evolution towards 
a peak structure in S(k) is observed near 2fcp. This peak 
is much enhanced at half filling, in agreement with BA. 

To summarize, when integral properties are addressed 
(local quantities, thermodynamics, total energy), the 
COM in the two pole approximation is accurate enough 
to yield a correct description of the system. The agree- 
ment with BA is indeed excellent in the case of half fill- 
ing and improves the results of other analytical methods 
which are more lengthy and more complex in many cases. 
The charge susceptibility is also well described as the 
charge excitations are dominated by the energy scale set 
by the opening of the Mott gap, and this is indeed caught 
by a two-pole approach. However, regarding the spin dy- 
namics the one-loop approximation can only grasp some 
general physics but it is too simple to investigate it prop- 
erly. We expect to receive better results whenever we will 
set up an approximation, for the two-particle propaga- 
tors, well beyond the one-loop one used here. Anyway, it 
is worth noting how this simple two-pole approximation 
is able to catch the spin and charge-dominated energy 
regions. 
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APPENDIX A: THE 2-POLE APPROXIMATION 

The doublet field (||) satisfies the Heisenberg equation 



4*(0 - - ( "? f (<) rn 4 rt ( il + A (i)1 ) ( A1 ) 

dt w ' \-(fJ,-U)r](t)+4tir(i) J y ' 
where 7r(i) is the composite field 

tt(z) = V c a (i) + c(i) [J a (i) c(i)] (A2) 

In the two-p ole a pproximation we linearize the equa- 
tion of motion (|Al|) as 



d 



(A3) 



where the energy matrix e{i, j) is defined by 

= E (I <{*(0, (A4) 
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In the two-pole approximation, by assuming transla- 
tional invariance, the thermal retarded Green's function 
S(k,u) = T [{11 [*(*) is given by 

S (k, u) = KyrHk) (A5) 



where s(k) = T [e(i, j)] is the energy matrix in momen- 
tum space and I(k) = T ({^(i),^^^')}) is the normal- 
ization matrix. 
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